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, Abstract 

] In the Minimal Supcrsymmctric Standard Model with bilinear _R-parity violation, only one 

• neutrino eigenstate acquires a mass at tree level, consequently experimental data on neutrinos 

, cannot be accommodated at tree level. We show that in the Next-to-Minimal extension, where 

' a gauge singlet superfield is added to primarily address the so-called /i-problem, it is possible to 

I generate two massive neutrino states at tree level. Hence, the global three-flavour neutrino data 
can be reproduced at tree level, without appealing to loop dynamics which is vulnerable to model- 

C^h! dependent uncertainties. We give analytical expressions for the neutrino mass eigenvalues and 

^ • present examples of realistic parameter choices. 

D 
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.p.', I Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) with biUnear interactions in the superpo- 
tential explicitly violating the i?-parity symmetry pP, a neutrino Majorana mass can be generated. 
Nevertheless, the rank 1 nature of the neutrino mass matrix suggests that only one eigenstate be- 
comes massive at tree level, whereas the neutrino oscillation data require at least two non-zero mass 
eigenvalues The bilinear i?-parity violating {flp ) soft masses induce one more non-zero eigen- 
value, but only at one-loop order. What happens if one considers the next-to-minimal version of the 
MSSM, called 'NMSSM' [3], in the presence of bilinear terms in the superpotential ? Here, the 
particle content is extended by one gauge singlet superfield. Our main result in this work is that two 
non-degenerate massive neutrino states now emerge at tree level. The upshot is that one can now 
reproduce the neutrino oscillation data with the superpotential parameters and gaugino masses just 
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from tree level physics. On the other hand, in the usual MSSM with bihnear ftp terms, the genera- 
tion of the second neutrino mass eigenvalue relies on the soft supersymmetry breaking scalar masses 
which in turn bring more uncertainties from the supersymmetry breaking mechanism; furthermore, 
uncertainties from loop dynamics creep in too. 

An increasingly important virtue of the NMSSM [3j (see ^ for phenomenological studies) is that it 
ameliorates the 'little hierarchy' problem originating from the requirement of large soft supersymmetry 
breaking scalar masses compared to the electroweak scale (in order to sufficiently push the lightest 
Higgs mass beyond the LEP limit). The NMSSM also provides a solution to the so-called /i-problem 
by arranging the vacuum expectation value (vev) of the gauge singlet scalar of the order of the 
supersymmetry breaking scale, so that the ^ parameter turns out to be at the electroweak scale. 

A recent paper 5^, in the context of NMSSM with ftp couplings, deals with the generation of neutrino 
masses where two eigenvalues arise at loop level. In another recent analysis 0, it has been shown 
that the NMSSM with bilinear ftp terms offers a possible mechanism of neutrino mass suppression, 
thus significantly reducing the hierarchy between /Xj and (defined below). Besides, an alternative 
extension of the MSSM explicitly breaking the i2-parity has been proposed in order to simultaneously 
address the /^-problem and the issue of small neutrino masses • 

In Section ^ we present the superpotential of the model we consider. In Section II I II we discuss the 
effective neutrino mass matrix. We present our numerical results in Section HVl Finally, we conclude 
in Section FVl 

II Superpotential 

The NMSSM superpotential contains two dimensionless couplings A and k in addition to the usual 
Yukawa couplings: 

Wnmssm = Y^^QiHuU^ + YfjQ.HdD^ + Yf^LiHaE^ + XSH^Hd + ^kS"" (1) 

where YJ^''^'^ are the Yukawa coupling constants (z, j, k are family indexes), and Qi, Li, Uf, D^, Ef, Hu, 
Hd, S respectively are the superfields for the quark doublets, lepton doublets, up-type anti-quarks, 
down-type anti-quarks, anti-leptons, up Higgs, down Higgs, extra singlet under the standard model 
gauge group. An effective /i term, given by X{s)HuHd, is generated via the vev of the scalar component 
s of the singlet superfield S. 

We now take note that in supersymmetric theories there is no deep underlying theoretical principle 
for the imposition of i2-parity as a symmetry However, there exist strong constraints on the 
fip couplings coming from various phenomenological considerations [51 ^] . Limits on neutrino masses 
and mixings have also been translated into tight upper bounds for couplings 

In the present paper, we consider a generic NMSSM superpotential containing both the bilinear and 
trilinear ^p terms: 

W = Wnmssm + IJ'iLiHu + XiSLiHu, (2) 

where /Xj (Aj) are the dimension-one (dimensionless) fip parameters. Actually, the contribution of 
trilinear term XiSLiHu was studied in Ref. [;5.. Admittedly, the most generic NMSSM superpotential 
also contains the other renormalizable trilinear fip interactions, namely, XijkLiLjE^, X'-f^LiQjD'^ and 
X'-jf^U[DjD'^, which are not relevant so long as we stick to tree level neutrino mass matrix. 
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Normally, in the NMSSM, only trilinear couplings with dimensionless parameters (like A and k) are 
kept in the superpotential, while dimensional parameters (like //) are generated from the vev (s). In 
this paper, the ftp NMSSM superpotential © is assumed to arise in either one of the following three 
possible scenarios: 

1. All possible renormalizable terms are included in the superpotential. Then both bilinear 
{fioHuHd, fiiLiHu) and trilinear {XSHdHu, XiSLiHu) terms are admitted. However, even if 
one may start with a term ^QHuHd, it can be rotated away by a redefinition of fields through a 
rotation on = (-ff^, Li) [a = 0, . . . , 3], since and have the same gauge quantum num- 
bers. There is no reason why this redefinition would remove also the Aq terms (Aq, = (A, Aj)), 
since the corresponding 4x4 rotation matrix depends on the parameters (the generic case 
is considered here, where fia and Aq, are not proportional). The coefficient fis of the 5^ term 
is assumed to be zero which can be considered as a possible natural value for a superpotential 
parameter. It should be noted that in the standard NMSSM with conserved i?-parity there is an 
accidental discrete symmetry whose spontaneous breaking causes the domain wall problem. 
In our version of the NMSSM, the fip bilinear LiHu term explicitly breaks that symmetry. 
In this scenario, we simply assume the existence of the dimensionful ^Uj terms, which as we will 
see later will be constrained from neutrino data. But we do not advance any theoretical reason 
as to why //j would be small. 

2. Our second scenario is based on the 't Hooft criteria of naturality: the parameters /Xj, /xq and fis 
are naturally small if the symmetry of the theory increases as these parameters are set to zero. 
For instance, one can assume that somehow a weak breaking (compared to the electroweak scale 
Qew) of some symmetry (like e.g. a U(l) symmetry forbidding the bilinear terms) generates the 
bilinear terms in the superpotential associated to fii, fio, fis ^ Qew- This small breaking would 
allow to address the ^-problem. Indeed, the main contribution to the dimension-one coefficient 
of HuHd here comes from A(s), as fiQ <^ fi = A(s) ~ Qew- The weak breaking of the symmetry 
is also responsible for the smallness of ftp couplings and neutrino masses, since //j <C Thus, in 
such a scenario, the treatment of the mii-naturalness (d la NMSSM) and of the neutrino masses 
(d la ftp ) are nicely connected via the weak breaking of a common symmetry. Admittedly, we 
do not provide any specific realisation of this weak breaking. We only hint at such a possibility 
that the bilinear fii, hq, fig couplings may arise from powers of some small spurion vev (<C Qew)- 

3. Finally, we propose a scenario where the trilinear Aj terms in superpotential are not present. 
This scenario relies on the Z^ symmetry, where the chiral superfields transform by exp(z27rg/3), 
with the following charge assignments: q = for C/^, D^, E'^; q = I for S, Hu, H^; and q = 2 
for Q, L. Such a symmetry allows all couplings in Eqs. and ^ except SLiHu- The other 
terms and HuHj, are also eliminated by this symmetry. A spontaneous breaking of this Z3 
symmetry admittedly creates the domain wall problem, as happens in the standard NMSSM 
with i?-parity. In this scenario, the /x term is created from the vev of S, but the /ij terms are 
present in the superpotential from the beginning only to be constrained by neutrino data. 

Ill Neutrino mass matrix 

Neutralino mass matrix: Within our framework, the neutralino mass terms read as, 

£™ = -^^'°^A^;^o^'° + h.c. (3) 
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in the basis ^'^^ = {B^, W^, h^, h^, s, i^i)^, where ^ {§) are the fermionic components of the super- 
fields ^ (S) and [i = 1, 2, 3] denote the neutrinos. In Eq. (jS)), the neutrahno mass matrix is given, 
in a generic basis (where (z>i) = Vi ^ 0, fj-i ^ and Aj 7^ 0), by 



-A^NMSSM 



03x3 



(4) 



where A^nmssm is the neutrahno mass matrix corresponding to the NMSSM. While writing the latter 
mass matrix, we assume Vi <C Vu,d (the exact expression of Mz being given by v'^ = ^^n + ^d + Z]i=i "^f — 
2cg M'^/g^ ~ (246/\/2 GeV)^). Also, we use s and c to stand for sine and cosine, respectively. 
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Above, Ml (M2) is the soft supersymmetry breaking mass of the bino (wino), tan/3 = Vu/vd = 
{h^)/{h^), and n = X{s). We assume for simplicity that A, k and the soft supersymmetry breaking 
parameters are all real. 

In Eq. ((IJ), is the part of the matrix mixing neutrinos and neutralinos: 

f -3^ 3^ ^i + Ai(s) Ait;„ \ 



V 



^ /Z3 + A3(S) X^Vu J 



(6) 



g and g' being the SU(2) and U(l) gauge couplings. 

Effective neutrino mass matrix: We restrict ourselves to the situation where Vi/v^^d 1 (as 
before), <C 1 and |Ai/A| <C 1 so that (i) no considerable modifications of the NMSSM scalar 

potential are induced by the additional bilinear and trilinear term in superpotential Q, (ii) the 
neutrino-neutralino mixing is suppressed, leading to sufficiently small neutrino masses as shown later, 
and (iii) the effective neutrino mass matrix can be written to a good approximation by the following 
see-saw type structure, 

mu = -i^^ -^NMSSM • C^) 

From Eqs. ©, (jEj) and Q, we deduce the analytical expression of the effective Majorana neutrino 
mass matrix: 

[^ii^jT + ifiiAj + fij A^)g + AiAjH] , (8) 



|-A4nmssm| 



where |A^nmssm| is the determinant of matrix ©, Aj = {s){Xi + A^) and: 



g 



\^vlMiM2 + X, 

3 

X + {Xvd - cos^ (3 XiVi) y, 



(9) 



i=l 
3 



n = X + 2008" l3{Xvd-Y,XiVi)y, 



(5) 



i=l 
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with, 



X = 2 cos^ (3 k{s)mI{cos^ Ow Ml ewM2), 
y = ^MUcos'^0wMi + sm'^9wM2). 

One should notice that the neutrino mass matrix (jSJ arises entirely at the tree level. 

The emergence of the specific mass matrix structure of Eq. ((HI at tree level is the primary result 
of our analysis. Accordingly to this particular structure, if the two (effective) quantities /Xj and Aj 
(characteristic of the NMSSM) take simultaneously non-vanishing values, then the mass matrix ceases 
to be of rank 1, even though the determinant is still zero. In this situation, we get two non-zero neutrino 
mass eigenvalues. It is worth comparing the situation with what happens in the MSSM with bilinear 
ftp violation. In the latter case, we get a similar kind of analytic structure of the mass matrix from 
the simultaneous consideration of the /ij as well as the corresponding soft Bi terms. While the f^ifij 
contributions originate at tree level, the i^iBj and BiBj contributions arise at one-loop order through 
Grossman-Haber diagrams [T7 which proceed through slepton-Higgs and neutrino-neutralino mixings 
(for a series of analysis in a three- flavour framework, see ^SI)- The Grossman-Haber loops would 
still contribute in our scenario, but now that we have two tree level masses, those loop-suppressed 
contributions are not so crucial for generating a viable neutrino mass spectrum. 

Neutrino mass eigenvalues at tree level: The eigenvalues of the effective neutrino mass matrix 
turn out to be {0,m~,m+} with. 



1 



m,. 



± < 



2 I NMSSM I 
/ 3 \ / 3 



yi=l / \i=l / \i=l / 



^4=1 



x + 



,1=1 



2\ 2 



(10) 



where T, g,n are defined in Eq. ©, and I = 4:{g^ - TH). 

Note that for either all /ij = or all Aj = 0, the eigenvalue m~ vanishes as expected since in this limit 
we recover the rank 1 form. An inspection of the form of Eq. (|1U() reveals that the coefficient of I can 
be written as Yli^jit^i^j ~ f^j^i)'^^ which indicates that the misalignment between /ij and Aj is crucial 
in creating a non- vanishing m~ . 

Therefore, the condition for generating two non-vanishing and non-degenerate eigenvalues is to ensure 
/ij 7^ and Aj 7^ simultaneously. In other words, to achieve two non-zero eigenvalues, fii has to be 
non-zero always, but we can go to a basis of Lq, fields where t^j = (then generally Aj 7^ 0), or to the 
other extreme where Aj = (but Vi 7^ 0), or any other basis in between basically maintaining Aj 7^ 0. 
On the contrary, if /ij = 0, only one neutrino eigenstate gets a mass different from zero, as was also 
found by the authors of Ref. But all the scenarios we considered in this paper, according to the 
above arguments, will yield two non-zero neutrino masses at tree level. We mention that within the 
first scenario, a rotation on the Lq fields has already been performed, whereas in the third one, no 
rotation is possible. 



In Figures (pQ) and ((21), we present the Feynman diagrams contributing to the Majorana neutrino mass 
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lajorana neutrino masses proportional to 



Figure 1: Tree level Feynman diagrams in the fip NMSSM generating 
hibj. The effective bilinear parameter hi stands for either /i^ or \i{s) . The mass-dimensional couplings appearing 
at the two vertices are of the type m — Mzt{6w) sin/3 [^(6*14^) — smOw for and t{6w) = — cosOw for W3] in 
(a) and m = —\vd in (b). A cross indicates either a mass insertion or a vev. The arrows show the flow of momentum 
for the associated propagators. 



Q. All these diagrams proceed through the tree level exchange of the neutralinos (gauginos, higgsinos 
and singlino). In these figures, we have considered the basis corresponding to = for simplicity. 



IV Numerical results 



Thus the present model predicts a hierarchical neutrino mass spectrum at tree level. This hierarchical 
pattern could be modified by the loop level contributions to neutrino masses. At least, the massless 
state acquires a mass once the loop contributions are considered. 

It is not our aim to give detailed numerical fits in this short paper. We would just like to numerically 
demonstrate that our scenario can reproduce the neutrino data from the tree level neutrino mass 
matrix, with a choice of NMSSM parameters that corroborate the //-naturalness. 

As an example, we take the following NMSSM parameters: 

A = 0.4, K = 0.3, /i = -200 GeV, tan/? = 30, Mi = 350 GeV, M2 = 500 GeV. (11) 

These parameters satisfy the general NMSSM constraints described below. These constraints are not 
expected to be significantly modified by the presence of fip interactions in superpotential as we 
work under the assumption ^\ <^ 1 and |Ai/A| <^ 1. 

The /i-naturalness forces one to restrict to (s) < 10 TeV, which translates into the condition [GeV] x 
10~^ < A. Furthermore, the absence of Landau singularities, for A, k, the top and bottom Yukawa 
coupling constants below the GUT energy scale, imposes |1] the typical bounds: A < 0.75, \k\ < 0.65 
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Figure 2: Tree level Feynman diagrams in the ftp NMSSM generating Majorana neutrino masses proportional to 
AiAjW^ and A,;u„6j. The effective ftp bilinear parameter bi, as in Figure 1, stands for either /i^ or Xi{s). The 
mass-dimensional coupling appearing at the right vertex in (b) is of the type m = —\vd- 



and 1.7 < tan /3 < 54. Finally, the LEP bound on the lightest chargino mass, namely > 103.5 GeV 
[Hj, translates into |^| > 100 GeV. 

Together with the values in Eq. ()11() . we take the following ftp effective couplings, 

fii = l- 10-^ GeV, fi2 = 2- 10-^ GeV, fis = 2 ■ 10"^ GeV 
Ai/(s) = 1.9 • 10-^ A2/(s) = 1.4 • 10-^ A3/(s) = 1.5 • 10"^ (12) 

This set of parameters yield the following three neutrino mass eigenvalues at tree level: 

m^i = 0, = 0.0095 eV, mj^g = 0.058 eV. (13) 

These values are in agreement with the three-flavour analyzes including results from solar, atmospheric, 
reactor and accelerator oscillation experiments which lead to {Aa level): 6.8 < Am'^i < 9.3 [10~^eV^] 
and 1.1 < Am|i < 3.7 [lO-^eV^] [2]. Besides, the neutrino mass eigenvalues in ()13() satisfy 
the bound extracted from WMAP and 2dFGRS galaxy survey (depending on cosmological priors): 
Yli=i "^fi ~ 0.7 eV [12]. Finally, these eigenvalues are perfectly compatible with the limits extracted 
from the tritium beta decay experiments (95% C.L.): m/j < 2.2 eV [Mainz] and mp < 2.5 eV [Troitsk] 
|16j . this effective mass being defined as m| = Yl^=i Weil'^ml. where Uei is the leptonic mixing matrix. 

Although we have chosen a particular set of input parameters for illustration, solutions exist over 
a wide range of parameter space. More realistic estimates can be obtained by switching on the soft 
scalar terms + h.c. plus the superpotential terms AjjfcLjLjE^^ and X'-jf^LiQjD'^. All these terms 

contribute to the neutrino mass matrix at one-loop order. A combined fit of all these parameters is 
beyond the scope of this paper. 



7 



V Conclusion 

In the NMSSM, which is a gauge singlet extension of the MSSM addressing the /i-naturalness, two non- 
vanishing neutrino mass eigenvalues can arise at tree level when the lepton number violating bilinear 
terms ^iLiH^ are present. One can then explain the neutrino oscillation data without essentially 
depending on the loop-generated masses which otherwise bring in more uncertainties from unknown 
dynamics. This result is in contrast with any other supersymmetric fip scenario, as those scenarios 
do not generate more than one massive neutrino eigenstate at tree level (except the scenario proposed 
in [7j where 3 right-handed neutrinos are added to the field content). 
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